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Abstract We consider in this paper a general two-sided jump-diffusion rislc model that 
allows for risky investments as well as for correlation between the two Brownian motions 
driving insurance risk and investment return. We first introduce the model and then 
find the integro-differential equations satisfied by the Gerber-Shiu functions as well as the 
expected discounted penalty functions at ruin caused by a claim or by oscillation; We also 
study the dividend problem for the threshold and barrier strategies, the moments and 
moment-generating function of the total discounted dividends until ruin are discussed. 
Some examples are given for special cases. 

AMS 2010 subject classifications: Primary 60J75; Secondary 60G51. 

Keywords: Paulsen-Gjessing's risk model; Stochastic return on investments; Integro- 
differential equation; Doleans-Dade exponential; Gerber-Shiu function; Dividends 



1 



1 Introduction 



The study of insurance risk models with stochastic return on investments has attracted 
a fair amount of attention in recent years, for example, Paulsen (1993) proposed the 
following general risk process Uf that allows for a stochastic rate of return on investments 
as well as a stochastic rate of inflation: 

The notation E{A) denotes the Doleans-Dade exponential of A given as the solution of 
the stochastic differential equation dE{A)t — E{A)t-dAt with E{A)o — 1, and Pt-,It and 
Rt are all semimartingales representing the surplus generating process, the inflation gen- 
erating process and the return on investment generating process, respectively. The initial 
values are Pq = u, Iq = and Rq — 0. He obtained an integro-differential equation and 
an analytical expression for ruin probability under certain conditions. Paulsen and Gjess- 
ing (1997a) simplified the model above by assuming that there is no inflation and both 
the surplus Pt and the return on investment Rt are independent classical risk processes 
perturbed by Brownian motions. Paulsen (1998a) considered a risk process Ut given by 

Ut^u + Pt+ [ Us-dRs, with Po = i?o = 0, (1.1) 

Jo 

where Pt and Rt arc independent Levy processes. With the above notation, the solution 
of (1.1) can be written as Ut = £{R)t{u + £ {R)j}dPs) . Cai and Xu (2006) considered a 
risk model that assumed the surplus of an insurer follows a jump-diffusion process and the 
insurer would invest its surplus in a risky asset, whose prices are modeled by a geometric 
Brownian motion. In the Discussion of the paper, Hailiang Yang extended the model of 
Cai and Xu to the case in which the surplus can be invested in both risky and risk-free 
assets. 

For some related discussions, among others, we refer the reader to Cai (2004), Yuen, 
Wang and Ng (2004), Cai and Yang (2005), Yuen and Wang (2005), Zhang and Yang 
(2005), Yuen, Wang and Wu (2006), Meng, Zhang and Wu (2007). For further references 
see two survey papers Paulsen (1998b) and Paulsen (2008). Some recent papers extended 
the model to renewal risk models with stochastic return, see e.g. Gao and Yin (2008) and 
Li (2012). 

Motivated by the previously mentioned papers, the aim of present paper is to generalize 
the model given in (1.1) by considering that Pt and Rt are general two-sided jump-diffusion 
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risk models that allow for risky investments as well as for correlation between the two 
Brownian motions driving insurance risk and investment return. The rest of the paper is 
organized as follows. In Section 2 we introduce the model. Integro-differential equations 
for the Gerber-Shiu functions are established in Section 3 and, in Section 4, we study 
the dividend payments under the threshold and barrier strategies. Finally, we give the 
concluding remarks. 



2 THE MODEL 



Assume that the surplus generating process Pt at time t is given by 

Np,t 

P^^u + pt + apWp,t-Y,Sp,i^ ^>0' (2-1) 

i=l 

where u is the initial surplus, p and ap are positive constants, {Wp^t}t>o is a standard 
Brownian motion independent of the homogeneous compound Poisson process X^S* ^P,i^ 
while {Sp^i} is a sequence of independent and identically distributed random variables. 
Unlike the model in Paulsen and Gjessing (1997a), we assume that Sp^i take values in 
(— oo, +oo). The upward jumps can be explained to be the random gains of the company, 
while the downward jumps are interpreted as the random loss of the company. Let Ap 
be the intensity of Poisson process Np^t, and Fp be the common distribution of Sp^i. We 
assume throughout this paper that EfSp^j] < oo and p — XpE[Sp^i\ > 0. 

Now, suppose that the insurer would invest its surplus a risky asset, whose price is 
assumed to follow the stochastic differential equation dS{t) — S{t—)dRt, where Rt is the 
return on investment: 

Rt = rt + aRWR,t + J2SR,i^ ^>0' (2-2) 

i=l 

where {WR^t}t>Q is another standard Brownian motion, independent of the homogeneous 
compound Poisson process Sn^i, while r and ap are positive constants. The intensity 

of Np^t is denoted by Xp, and the distribution function of the jump 5*^ by Fp. Unlike 
the model in Paulsen and Gjessing (1997a), Wp^t is correlated with Wp^t and Wp^t can be 
written as Wp^t = pWp^t + a/1 — p^Wp^, where p G [—1, 1] is constant, Wp^ is a standard 
Brownian motion independent of Wp^f When = 1, there would only be one source of 
randomness in the model. 
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We define the risk process Ut as the total assets of the company at time t under this 
investment strategy, then Ut is the solution of the stochastic differential equation 

dUt ^ dPt + Ut-dRt, t>0. (2.3) 

By using Theorem 1 in Jaschke (2003) it is not hard to see that the solution of (2.3) is 
given by 

Ut = Si{R)t {^ + £i{R):hdPs - P<^p<^rJ^ ^i(^)rids^ , (2.4) 

where 

Si{R)t = exp I (^r - t + aRWn^t^ f[{l + Sr,^). 

Because the quadratic variational processes of 

(TpWp,t+ [ aRUs-dWR,s^ [ {ap + paRUs-)dWp,s+ [ apy/l - p''Us-dW°^, 
Jo Jo Jo 

and ^ 

^ ^Ji<Jp + PCJrUs-Y + 4(1 - P'WL^B, 

are same, where {B}t>o is a standard Brownian motion independent of the compound 
Poisson processes involved, by Ikeda and Watanabe (1981, Theorem 7.2, p. 85), they have 
the same distributions. Thus, in distribution, (2.3) can be written as 

Ut^u + f {p + rUs)ds 
Jo 

- I Us-d\Y,SRA . (2.5) 



Using Ito's formula for semimartingale, one finds that the infinitesimal generator C of 

{Ut}t>o is given by 

^9{y) = ^(4 + 2pap(7Ry + 4?/')/(?/) + (P + ry)g'{y) 

/oo 
[9{y -z)- g{y)]dFp{z) 
-oo 

/oo 
^[g{y + yz)-g{y)]dFR{z). (2.6) 

Reiricirk 2.1. If Fp{0) — and p — 0, then (2.4) reduces (2.4) of Paulsen and Gjessing 
(1997a) and (2.6) becomes (2.6) of Paulsen and Gjessing (1997a). 
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3 The Gerber-Shiu functions 



In this section, we consider the Gerber-Shiu expected discounted penalty function for the 
risk process (2.4). The time of ruin of (2.4) is defined as T — mi{t > : C/j < 0} with 
T = oo if ?7t > for alH > 0. The ruin probabihty with an initial surplus m > is defined 
as 

iIj{u) = P(r < oo\Uo = u). 

Note that ruin may be caused by a claim or by oscillation. Denote the ruin probabilities 
in the two cases by 

ips{u) = P{Ut <0,T < oo\Uo = u), ^d{u) = P(C/t = 0,7 < oo|;7o = u). 

Obviously we have 

'il){u) = tlJsiu) + i)d{u), u>0. 
Moreover, when ap ^ 0, it follows from the oscillating nature of the process Ut that 

V'(O) - V'd(O) = 1 and V.(0) = 0. 

Let w = w{xi, X2) be a nonnegative bounded measurable function on [0, 00) x [0, 00). As 
in Gerber-Shiu (1998), the Gerber-Shiu expected discounted penalty function is defined 

by 

(j){u) = E[e-^^w{UT-, \Ut\)I{T < oo)\Uo = u], (3.1) 

where 5 > is a constant and I{A) is the indicator function of event A. Sometimes we 
write Pu for the probability law of U when Uq — u and E„ for the expectation with respect 
to P„. 

Following Wang and Wu (2008), we decompose the Gerber-Shiu function 0(m) in (3.2) 
correspondingly into the following two parts: 

4>s{u) = E[e-'^w{UT-, \Ut\)I{Ut < 0,T < oo)|t/o = u], (3.2) 
(t)^{u) = 'u;(0, 0)E[e-'^^/(C/r = 0, T < 00) |C/o = «]. (3.3) 
Obviously, 'ip{u), i)s{u) and '4^d{y) are special cases of 0(w), 0s(w) and 4)d{u), respectively. 

For simplicity, we define the operator 

Qh{u) = -(c7"p + 2papaRU + a\u^)h" {u) + (p + ru)h'{u) 

/u roo 
h{u - z)^Fp{z) + Xr h{u + uz)^Fr{z). (3.4) 
-00 J -I 
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Theorem 3.1. Assume that (f){u),(f)s{u),(f)d{u) are twice continuously differentiahle on 
[0,00), where the derivative at u = means the right-hand derivative. If aj, > 0, p — 
pcpCR ~ ApE[S'p,j] > 0, 1) = and r — \(t\ > 0. Then, 

(i) (j){u) satisfies the integro- differential equation 

POD 

{S + Xp + Xr)(P{u) = g(p{u) + Xp w{u,z-u)dFp{z), u>Q, (3.5) 

J u 

with boundary conditions 

0(0+) = w{0, 0), lim (f){u) = 0. (3.6) 

w— >-oo 

(ii) 0s (u) satisfies the integro- differential equation 

POO 

{S + Xp + XR)(psiu) = g(psiu) + Xp w{u,z-u)dFp{z), u>Q, (3.7) 

J u 

with boundary conditions 

0,(0+) = 0, lim (j),{u) = 0. (3.8) 

(iii) (pdi'U') satisfies the integro- differential equation 

(5 + Ap + Aij)0d(w) =e;0rf(ti), u>0, (3.9) 

with boundary conditions 

0d(O+) = w{0, 0), lim = 0. (3.10) 

Proof. We can prove the theorem following a similar argument as in Yin and Wang 
(2010) by using Ito's formula. In the following proof, however, we using a more intuitive 
infinitesimal argument as in Cai and Yang (2005) and Cai and Xu (2006), where the ruin 
probabilities have been studied. The main difference from theirs is that our model has 
two Poisson processes and has two-sided jumps. 

(i) Let 

Yt = ue^' + pe^' I e-^^ds + ape^' [ e-^^dWp,„ (3.11) 
Jo Jo 



where 



At = { r - ^4 ) t + apWR^t. 



Consider the risk process Ut, defined by (2.4), in an infinitesimal time interval (0, t]. Since 
both Np^t and Np^t are Poisson processes, there are five possible cases. 
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(i) . Both Np^t and A^^ ^ have no jumps in (0,t] (the probabihty that this case occurs 
is e-^^*e-^«*). Thus Ut = Yt. 

(ii) . There is no jump of Np^t in (0,t] and there is exactly one jump of Np^t in (0,t] 
(the probabihty that this case occurs is e~'^^*Apte~^^*), with claim amount z, and 

(a) z <Yt, i.e. ruin does not occur and, thus Ut = Yt — z. 

(b) z > Yt, i.e. ruin occurs due to the claim, or 

(c) z — Yt, i.e. ruin occurs due to oscillation (the probability that this case occurs is 



(iii) There is no jump of Np^t in (0, t] and there is exactly one jump of Nji^t in (0, t] 
(the probability that this case occurs is e'^^^Xpte"^^*), and thus Ut — {1 + SR^i)Yt. 

(iv) Both Np^t a-nd Np^t have one jump in (0, t] (the probabihty that this case occurs 



(v) Np^t and/or Np^t has more than one jumps in (0, t] (the probability that this case 
occurs is o{t)). 

By considering the five possible cases above and noticing that in case (ii)(b), (t){Yt—z) — 
w{Yt, z — Yt), we have 



zero) . 



is o{t)). 



4>{u) 



^-St^-Xpt^-Xnt-^^^^Yt) 

+e-'^*(Apte-^^*)e-^«*E, 



'u 





+e-^*(Apte-^^*)e-^«*E, 



'u 



w{Yt,z-Yt)dFp{z) 




'u 




(3.12) 
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Note that e-^*e--^-p*e-^«* = I - [S + Xp + Ar)^ + o{t), we have 
= E^<f>(Yt) - (l)(u) - (5 + Xp + XR)tE^<f>(Yt) 

f-Yt 

J —oo 
poo 

+e-**(Apie-^^*)e-^«*E„ / w{Yt, z - Yt)dFp{z) 

JYt 

/oo 

(j){Yt{l + z))dFR{z)+o{t). (3.13) 

By Ito's formula, we have 

lim ^ 1 (^2 ^ 2papaRU + alu^)(t>"{u) + (p + rM),^'(M). (3.14) 

Therefore, by dividing t on both sides of (3.13), letting t ^ 0, and using (3.14), we get 
(3.5). Let 

Jo Jo 
and Zoo — l™t->oo ^t- Since p — papaR — XpE[Sp^i] > 0, it follows from Theorem 3.1 in 

Paulsen (1993) that Zt is a submartingale. In addition, the conditions ap > 0, Fr{—1) = 

and r — |(7p > imply that Z^c exists and is finite with probability one. Thus, from (2.4) 

we find that Ut drift to +oo with probability one. Consequently, '0(+oo) = 0. The 

boundary condition lim„^oo 0(^) = follows from (p{u) < Mipi^u), where M is a upper 

bound of w{xi,X2); The boundary condition 0(0+) = w{0,0) follows from the oscillating 

nature of the sample paths of Ut- The results (ii) and (iii) can be proved by the same 

arguments as (i). 

Remark 3.1. Let us compare our results with known results. 

(1) . Letting p — 0,5 — 0, Fp{0) — and w{x,y) = 1 in (3.5) we get the result Theorem 
2.1 (i) in Paulsen and Gjessing (1997a); Letting p — 0, Fp{0) — and w{x,y) = 1 in 
(3.5) we get the result Theorem 2.1 (ii) in Paulsen and Gjessing (1997a). 

(2) . Letting p = 0, (Tij = 0, Ap = 0, Fp(0) = and w{x,y) = 1 in (3.5), (3.7) and (3.9), 
we get the result (3.1), (3.4) and (3.9) in Cai and Yang (2005), respectively. 

Example 3.1. Under the assumptions of Theorem 3.1, assume that 5 = Xp = Xr = 
and w{x, y) = 1, then for any u > 0, ip{u) and i/jdiu) satisfy the same differential equation 

i(o-^ + 2papaRU + a%u'^)h"{u) + {p + ru)h'{u) = 

and the following boundary conditions 

/i(0) = 1, /i(+oo) = 0, 
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where h{u) = iIj{u) or ipd{i^)- 



If IpI < 1, the solution was found by Haihang Yang; See Cai and Xu (2006, p. 130). If 
IpI = 1, the solution is given by h{u) = 1 — where 



2r 



Example 3.2. Under the assumptions of Theorem 3.1, assume that \p — Xr = 0, w{x, y) = 
1 and 5 > r, then for any w > 0, (j){u) and 4)d{u) satisfy the same differential equation 

^((Tp + 2papaRU + a]iU^)g"{u) + (p + ru)g\u) = 6g{u) (3.15) 

and the following boundary conditions 

5(0) = 1, 5(+oo) = 0, 

where g{u) = (f){u) or (pdiu). 

A change of variables x = u + —p and h(x) = giu) brings the equation (3.15) into 
the form 

\ [4(1 - p") + alx^] h"{x) + n{x)h\x) = 6h{x), (3.16) 



where 



/ N rpop 
K[x) — p h rx. 



When < 1, (3.16) has the same form as (Al) in Paulsen and Gjessing (1997a), using 
Theorem A.l. in Paulsen and Gjessing (1997a) we get 

h{x) = CiD{x, a + l) + C2E{x, a + 1), 

where 

/Jctan((aH/ap)x)(cost)^-^(apsint - a rx COS t)'' exp {-^t} dt 



D{x,X) 



E{x,\) 



/^|*"°^^"«/"^^"^(cost)^-^(a^xcosi - apsint)^exp {-^t} dt 



Here 
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Because E(x, a + 1) — >■ +00 as a; — >■ +00 and h{+oo) — g{+oo) — 0, so that C2 = 0. In 
addition, using boundary condition ^'(0) = 1, we get 

Thus, 

g{u) = h{x) = 

This result is also obtained by Paulsen and Gjessing (1997a) in the case where p = 0. 



4 Total discounted dividends 



4.1 Threshold strategy 

In this subsection, we consider the threshold strategy for dividend payments. More specif- 
ically we assume that the company pays dividends according to the following strategy 
governed by parameters 6 > and // > 0. Whenever the modified surplus is below the 
threshold level 6, no dividends are paid. However, when the surplus is above this threshold 
level, dividends are paid at a constant rate p. Once the surplus is negative, the company 
is ruined and the process stops. We assume that the risk process U without dividends 
follows (2.3). We define the modified risk process C4 = {Ub{t) : t > 0} in which Uij{t) is 
the solution of stochastic differential equation 

dUbit) = dPt + Ub{t-)dRt - dDbit), t > 0. 

where Di,{t) — ji J^I{Ui,{s) > b)ds. Let Di{b) denote the present value of all dividends 
until time of ruin Ti, 

Jo 

where Ti = inf {t > : Ub{t) < 0} with Ti = 00 if Ub{t) > for all t > 0. Here 5 > is the 
discount factor. Denote by V{u; b) the expected discounted value of dividend payments, 
that is, 

V{u) := V{u;b) = E[Di(6)|[/6(0) ^u] = E„[L'i(6)]. 
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Let 

Mi{u,y] b) = E[e^^i(*)|[/6(0) = u] = E„[e^^i(*)] 

denote the moment-generating function of Di{b). If 5 > 0, then < Di{b) < |, and thus 
M{u, y; b) exists for all finite y. 

Theorem 4.1. Assume that V{u) is twice continuously differentiahle on (0, &) U (6, cxo). 
Then for < u < b, V{u) satisfies the following integro- differential equation: 

{5 + \p + XR)V{u)^gV{u), (4.1) 

and for u > b, V{u) satisfies the following integro- differential equation: 

{S + Xp + Xr)V{u) = gV{u) - fiV'iu) + fi, (4.2) 

where Q is defined by (3.4)- 

Proof. Let V^(h) be twice continuously differentiahle and equals to V{u) on (— oo, 6 — 
^] U [6-1- ^, oo). Applying Ito's formula for semimartingales to deduce that for t e [0, Ti) 

e-^'V^mt)) = V^mO)) + f e-''{C^ - 5)Vmms))ds + M^, 

Jo 

where M'^ is a local martingale and is defined as C^g{y) = —iJ.I{y > b)g'{y) + Cg{y), 
where C is defined by (2.6). It follows that for any appropriate localization sequence of 
stopping times {Tn,n > 1} we have 

r ptATiATn 

E^[e-'(^^'^^^-r.)Y^(Ub{tAT^Arn))]^Vm{u;b) + E^ / e-'^L^ - 5)Vm{Ub{s))ds . 

Jo 

(4.3) 

Letting n, m t oo and 1 1 iii (4-4) and note that V{Ub{Ti)) = 0, we find that 



V(u) = E„ 



r-Li 

/ ne~^'I{Ubis) > b)ds 
Jo 



if and only if C^V{u) — SV{u) — —^I{u > b). From which we get (4.1) and (4.2). This 
ends the proof of Theorem 4.1. 

Remark 4.1. It can be verified that V{u) = on it < 0, V{Q) — if ap > and 

limm.oo V^('u) = f; V satisfy the continuity condition V{b—) = V{b+) = V{b). Moreover, 
ifap = aR = thenpV'ib-) = {p - ^i)V'{b+) + n, and if a p > 0, then V'{b-) =V'{b+). 

Remctrk 4.2. The above result was obtained by Gerber and Shiu (2006) for the compound 
Poisson model, Wan (2007) for the compound Poisson model perturbed by diffusion and 
Ng (2009) for the dual of the compound Poisson model. 
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Example 4.1. Assume that Xp = Xr = 0. Then V{u) solves the following different 
equations 

^(cTp + 2papaRU + a\u^)V"{u) + (p + ru)V'{u) = 5V{u), < u < b, 

^((Tp + 2papaRU + a\u^)V"{u) + {p- ii + ru)V'{u) + = 5V{u), u > b, 
with the boundary conditions 

V{0) = 0, lim V{u) = ^, V{b-) = T/'(6-) = V\b+). 



Similar to Example 3.2, the solution is given by 



CsD{u + ^p, a + 1) + C4E{u + ^p, a + 1), ii u < b, 
CMu + £|p, a + 1) + C,E,{u + ^p, a + 1) + f , if M > b, 



Di(x, X) = 



where D,E, a and /3 are defined in Example 3.2 and 

/Jctan((.W^pW(^Q^^)^"'(^^^^"^ - ^il^COSt)^exp {-f^t} dt 

ap^+^iaRy+>^ 

/:f ""^^"«/""^"^(cosi)/'-^((7pa;cost - apsint)^exp j-M^j dt 
^^^^'^^ = ^ api-H/^(a«)iH-A ^^^^^ ■ 

The constants C3 — Cg can be determined by the boundary conditions above and they are 
given by Ce = 0, 

^ _ fD,{b+ff,a)E{ff,a + l) 



^3 



Q,E{^,a+l)-Q,D{!^,a+iy 



^ fi^(^.n + l)D,(/; + ^.n) 

' Q,E{^,a + l)-Q,D{f^,a + iy 



Here 



Q^^D{b+P^,a+l)D,{b+^,a)-D{b+^,a)D,{b+P^,a+l), 

C^R CTr (Tr Gr 

Q, ^ D,{b + a)E{b + a + 1) + + « + + a). 

Theorem 4.2. Assume that Mi{u, y; b) is twice continuously differentiable in u on (0, 6)U 
(6,00) and once in y >0. Then Mi satisfies the following integro-differential equations 

AMi{u, y; b)-6y ^^'^^^ ^' -{Xr+Xp)Mi{u, y; b)+Xp{l-Fp{u)) = 0,0 <u<b, (4.4) 
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and 

AA/r ( u\ dMi{u,y;b) dMi{u,y;b) 
AMi{u,y;b) - n — 6y — h iiyMi{u,y;b) 

- {XR + Xp)M^{u,y;b) + Xp{l-Fp{u)) = 0, u>b, (4.5) 

where 

A nit h\ 1^ 2^0 ^2 2^9'^M{u,y;b) dM{u,y;b) 
AM{u, y; b) = 2^^^ 2papaRU + dj^u ) — h (p + ru) — 

/u 
M{u-z,y;b)dFp{z) 
■00 

/oo 
^ M{u + uz,y;b)dFR{z). (4.6) 

In addition, Mi {u, y; b) satisfies 

Mi(0,t/;6) = 1, (4.7) 
\im Mi(u,y;b)^ey^'/^. (4.8) 

Proof. When < u < b, consider the infinitesimal time interval from to t. By the 
Markov property of the process Ut, we have 

Ml («,!/; 6) = E„[e^/"^-"'^^W]+o(i) 

= E„[e^/o"-*^-^^*^^''^^(*+^)] + o(t) 

= E„[e^/"^-*'*+='^^Wo^,]+o(t) 

= E„(E,;Je^-"/o"-"<^^(^)]) + o(t) 

= E„[Mi(C/t,ye-^*;6)] + o(t), (4.9) 

where 9t is the shift operator. We refer to Kallenberg (2006) for more details on the 
Markov property and the shift operator. By the law of double expectation, we have 

E^[Mi{Ut,ye-''-b)] = (1 - Apt)(l - A«t)E,[Mi(y/, ye-^*; 6)] 

+Xpt{l - XRt)E^ / Mi(y/ - z,ye-''; b)dFp{z) 
J —00 

+Xpt{l - ARt)E,(l - Fp(F/)) 

/oo 
Mi{Y,\l + z),ye-'';b)dFR{z) 

+o(t), (4.10) 
where ^ ^ 
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Here 



By Ito's formula and note that 



Yl^u + I {p + rYs)ds 
Jo 

we have 

ou oy 
+-(4 + 2p<jp<7«F/_ + aliXlf) ^ 



Thus 



_,,E.]^-^^^«^p^d. (4.11) 

Substituting (4.11) into (4.10) and then dividing both sides of (4.10) by t, letting t ^ 
and rearranging, we obtain (4.4). 

Similarly, when u > b we have 

M,{u,y;b) = ey'''E^[M^{Ut,ye-''■,b)]+o{t), (4.12) 

from which we obtain 

Mi{u,y;b) = {l-Xpt){l-XRt)ey'''E4Mi{Y,^ye-'';b)] 



+Xpt{l 


- XRt)ey'''E^ J 


M^{Y^ - z,ye-'';b)dFp{z) 

-oo 


+Xpt{l 


- XRt)ey^'E,(l 


- Fp{y')) 


+ XRt{l 


-Xpt)ey>''E^J 


•oo 

M,{Y,\l + z),ye-'';b)dFR{z) 


+o{t), 




(4.13) 



where 



JO JO 
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Using the same argument as for (4.4) we get (4.5). The condition (4.7) is obvious, and 
(4.8) follows from limu^.(x, ^i(^) = f • This ends the proof of Theorem 4.2. 

Set 

oo ^ 

M,{u, y;b)^l + J2 yVkiu; b), (4.14) 
k=l 

where 

Vk(u) = Vkiu; b) = E[D,(b)''\Ub(0) = u], (4.15) 

is the kth moment of Di{b). Substitution of (4.14) into (4.4) and (4.5) and comparing 
the coefficients of y'^ yields the following integro-differential equations 

gVk{u) = {Xr + \p + kS)Vk{u), 0<u<b, (4.16) 

and 

gVkiu) - ^^^^ + ki^Vk-iiu) = (Xr + + k5)Vkiu), u > b, (4.17) 

where Q is defined by (3.4). They generalize (4.1) and (4.2), which are for k — 1. The 
boundary conditions are Vk{0;b) = and limu^^Vk{u;b) = (^)^. 



4.2 Barrier strategy 

It is assumed that dividends are paid according to a barrier strategy ^b. Such a strategy 
has a level of the barrier b > 0, when the surplus exceeds the barrier, the excess is paid 
out immediately as the dividend. When the surplus is below b, nothing is done. Let 
be aggregated dividends up to time t by insurance company whose risk process is modeled 
by (2.3). The controlled risk process when taking into account of the dividend strategy 
is — {Uf : t > 0}, where is the solution of stochastic differential equation 

dUf = dPt + U^dRt - dDl t > 0. 

Denote by Vi{u; b) the dividend- value function if barrier strategy is applied, that is, 

t/i(M;6) = E„[D2(6)], 

where D2{b) = e~^*dD^. Here J > is the force of interest for valuation and T2 ~ 
mi{t >Q:U^ <Q}. Let 

M2(«,y;6)=E4e^^^('')] 
denote the moment-generating function of D2{b). M{u,y,b) exists for all finite y. 
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Theorem 4.3. Assume that M2{u,y; h) is twice continuously differentiahle in u on (0,6) 
and once in y >0. Then M2 satisfies the following integro- differential equation 

AM^iu, y- h) - Sy^M^^^iyiR _ (A^ + Xp)M2{u, y; b) + Ap(l - Fp{u)) = 0, < m < 6, 

oy 

(4.18) 

where A is defined by (4-6). In addition, M2{u,y;b) satisfies 

M2(0,y;6) = l, (4.19) 
dM2{u,y-b)^ 



du 



yM2{b,y;b). (4.20) 



Proof. The proof of (4.18) is same as the proof of (4.4). The condition (4.19) is 
obvious. To prove (4.20), we first consider the special case in which ap — ap — 0. Letting 
It t 6 in (4.18) gives 

dM2{u,y;b) dM2{u,y; b) 
[P + rb) — \u=b - oy — \u=b - [Xr + Ap)M2(6, y; b) 

+Ap / M2{b - z, y- b)dFp{z) + Xp{l Fp{b)) 
J —00 

/oo 
M2{b{l + z),y;b)dFR{z)^0. (4.21) 

Similarly, for u — b we have 

M2{b, y; b) = e-^^*e-^«*e^(^'+"'')*M2(6, ye'^'; b) 

,b 

+Apte-^^*e-^«*e^(^+''''^* / M2{b - z, ye'^'; b)dFp{z) 

J —00 

+Apte-^^*e-^«*e^(P+^'')*(l - Fp{b)) 

/oo 
M2(6(l + z),ye-^*;6))dFp(z) 

+o{t). (4.22) 

This, together with the following Taylor's expansion 

-st.u._n.^u .,.5M2(6,y;6) 



gives 



M2(6, ye-''; b) = M2(6, y; b) - Syt ^^^^ + o{t). 

oy 



^^ dM2{u,y;b) ^^^^ _ ^^^ + Xp-yip + rb))M2ib,y;b) 



dy 



+Ap f M2{b-z,y;b)dFp{z)+Xp{l-Fp{b)) 
J —00 

/oo 
^ M2(6(l + ;2),y;6)dFp(z) =0. (4.23) 
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Comparing (4.23) with (4.21) we obtain 

dM2{u,y; h) . 



yM2{b,y;b). 



du 

If CTp > 0, (7^ = or cTp > 0, (7^ > 0, the process can be viewed as the hmit of a family of 
same kind processes without Brownian motions, and in this way (4.20) can be obtained 
as hmiting results. This ends the proof of Theorem 4.3. 

Set 

fe=i 

where 

is the kih moment of D2{b). Substitution of (4.24) into (4.18) and comparing the coeffi- 
cients of y'' yields the following integro-differential equation 

GVk(u) = (Xr + Xp + k5)Vk(u), 0<u<b, 

where Q is defined by (3.4). The boundary conditions are Vfc(0; b) = and 

dVk{u;b)^ 



du 



\u=b = kVk-i{b; b). 



Example 4.2. Assume that Xp — Xr — 0. liu <b, then Vi{u; b) solves the following 
differential equation 

- (4 + 2papaRU + alu^) ' + (p + ru) ' = 5V^ (u; b) (4.25) 

with 

V',(0;6)=0, ^^U.= l. (4.26) 
When < 1, the solution of (4.25) is given by 

Vi{u- b) = CrD lu + ^p, a+lj + C^E lu + ^p, a + 1 j , 



where D, E and a are defined in Example 3.2. The constants C-j and Cg can be determined 

^D{y,a) = -aD{y,a-l) and 



by conditions (4.26). Using that ^D{y,a) = —aD{y,a — l) and ■^E{y,a) = aE{y,a — l) 



we obtain 

{a + l)A{b+^p,ay 
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where 



(a + l)A{b + ^p,ay 



A(b+—p,a]^ E (b+—p,a] D (—p,a + l 
cfr J \ cfr J \<7r 

+D (b + — p, a] E i—p, a + 1 
\ ctr J Kctr 



It follows that 

B(u + ^p,a 



{a + l)A{b+^p,ay 

where 

B(u + —p,a]= D(—p.a + l]E(u+—p,a + l 
\ o-R J \aR J \ (Tr 

-E ( — p, a + 1 ) D ( w + — p, a + 1 

In particular, when p = we recover the result of Example 2.2 in Paulsen and Gjessing 
(1997b). 



5 Concluding remarks 



In this paper, a generalized Paulsen-Gjessing's risk model is examined, some rather gen- 
eral integro-differential equations satisfied by the Gcrber-Shiu functions, the expected 
discounted dividends up to ruin and the moment generating functions of the discounted 
dividends arc presented, respectively. Generally speaking, it is difficult to find the ana- 
lytical solutions except for some specials. A numerical method called the block-by-block 
has been used by Paulsen et al. (2005) to find the probability of ultimate ruin in the clas- 
sical risk model with stochastic return on investments. The solutions, either analytical 
or numerical, of the integro-differential equations in this paper are not only interesting 
but also valuable in practice. Other research problems such as the optimality results for 
dividend and investment need also to be studied. 
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